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ABSTRACT 

Analyses  and  interpretation  of  surface  layer  and  synoptic -scale  data 
obtained  in  the  Northeast  Atlantic  were  performed  to  obtain  descriptions 
of  the  evaporation  duct  and  associated  atmospheric  and  oceanic  synoptic 
features.  The  surface  layer  data  were  quite  unique  because  they  were  ob- 
tained from  high  quality  measurements  from  ships  spatially  separated  in  a 
fixed  array.  Magnitudes  and  horizontal  homogeneity  of  duct  heights  were 
compared  to  air-mass  trajectories  and  weather  patterns.  The  mean  duct 
height  was  4- 5m  and  higher  values  (8-13m)  occurred  with  North  and  North- 
west trajectories.  During  times  when  the  heights  were  8- 13m,  horizontal 
homogeneity  occurred  35%  of  the  time  with  a  maximum  duration  of  12  hours. 
This  was  established  on  the  basis  of  a  two  hundred  kilometer  separation 
between  locations  of  duct  height  estimates. 

Significant  features  of  the  duct  height  and  its  temporal  and  spatial 
variations  were  related  to  synoptic  scale  descriptions.  Satellite  imagery 
used  in  conjunction  with  point  observations  appears  to  provide  the  most 
useful  information  in  describing  the  intensity  and  areal  distribution  of 
the  evaporation  duct. 
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I.   INTRODUCTION 

Over  the  surface  of  the  ocean  there  exists  a  phenomenon  described 
as  the  evaporation  duct.  The  military  significance  of  its  existence 
has  been  known  for  many  years.  Attempts  have  been  made  to  quantify 
the  evaporation  duct  height  and  determine  its  effect  on  electro- 
magnetic propagation.  Assessment  of  evaporation  duct  characteristics 
in  tactical  military  applications  should  only  require  routine  ship- 
board measurements  and  simple  calculations;  hence,  a  great  deal  of 
study  has  been  performed  to  parameterize  the  duct  on  the  basis  of 
simple  observational  measurements  of  surface  and  near-surface  para- 
meters. Also,  as  the  use  of  low-flying  supersonic  aircraft  and 
missiles  become  increasing  threats  to  Naval  surface  combatants,  timely 
and  precise  descriptions  of  the  nature  of  the  evaporation  duct  will  be 
required  to  permit  its  exploitation. 

It  has  been  expected  that  in  an  area  in  the  high  latitudes ,  well 
away  from  large  land-mass  influences ,  and  where  frequent  migratory 
weather  systems  disturb   the  nature  of  the  atmosphere  near  the  surface, 
the  height  of  the  evaporation  duct  is  likely  to  be  small  and  its 
occurrence  highly  localized.  Past  research  has  been  primarily  limited 
to  regions  where  ducts  are  known  to  be  both  strong  (high)  and  per- 
sistent (e.g.  tropics,  southeast  portion  of  subtropical  high,  coastal 
regions  of  warm  offshore  winds)  with  relatively  little  examination  of 
higher  latitude,  open  oceanic  regions.  This  study  will  examine  the 
occurrence  and  nature  of  the  evaporation  duct  over  the  North  Atlantic 


Ocean  subject  to  highly  varying  atmospheric  conditions.  This 
examination  will  focus  on  determining  significant  ducts  and  their 
1)  frequency  of  occurrence  and  duration,  2)  temporal  and  spatial 
nature,  and  3)  sensitivity  to  atmospheric  conditions. 

The  results  reported  here  were  made  utilizing  a  unique  set  of  data 
acquired  from  the  Joint  Air-Sea  Interaction  Project  (JASIN)  conducted 
in  the  Northeast  Atlantic  during  the  late  summer  of  1978.  The  results 
were  analyzed  and  comparisons  were  made  between  several  ships ,  syn- 
optic analyses  and  high  resolution  satellite  imagery.  The  resultant 
examination  yielded  a  description  of  an  evaporation  duct  character- 
istic of  the  region  as  it  was  affected  by  frequently  changing  weather 
patterns . 

This  study  was  on  data  collected  from  four  research  vessels 
engaged  in  the  JASIN  experiment  from  mid- July  to  early  September,  1978 
under  open  ocean  conditions.  The  quantity  of  data  and  the  large  area 
of  study  permitted  an  examination  to  determine  the  nature  of  the  duct 
and  the  atmospheric  disturbances  most  likely  to  influence  it.  An 
important  aspect  of  the  experiment  was  the  positioning  of  four  ships  in 
a  triangle  (180  x  180  x  200  km)  during  periods  of  intensive  meteorologi- 
cal measurements .  This  constant  spatial  separation  permitted  accurate 
assessments  of  spatial  as  well  as  temporal  variations  over  a  large 
horizontal  region  as  weather  systems  moved  through  the  area.  The  data 
were  used  with  expressions  presented  by  Fairall  et  al  [1978]  to  compute 
the  associated  oceanic  evaporation  duct  height. 


II.  BACKGROUND 

It  has  been  long  recognized  that  distributions  of  temperature  and 
humidity  in  the  lower  atmosphere  are  responsible  for  anomalous  electro- 
magnetic (EM)  wave  propagation   [  Katzin,  1947  and  Kerr,  1951  ]   result- 
ing in  either  enhancement  or  degradation  of  Naval  sensors  and  weapons 
systems.  In  particular,  propagation  between  terminals  reasonably  close 
to  the  sea  surface  has  been  observed  to  exhibit  significant  variation 
of  signal  strength  at  ranges  far  in  excess  of  the  radar  horizon. 
(Fig.  1)  Such  anomalous  propagation  effects  were  first  seriously 
studied  during  and  immediately  following  World  War  II.  In  recent  years, 
advances  in  mixed- layer  theory  have  led  to  numerous  computational 
approaches  for  real-time  assessment  of  atmospheric  refractive  effects. 

Often  relatively  cool  moist  marine  air  extends  vertically  from 
the  ocean  surface  to  altitudes  of  a  few  hundred  meters.  Above  these 
heights  the  air  becomes  much  warmer  and  drier  for  a  variety  of  reasons. 
Therefore,  a  transition  layer  exists  in  which  the  air  temperature  in- 
creases and  moisture  rapidly  decreases  with  increasing  height.  The 
resulting  gradients  can  cause  the  EM  refractive  index  to  decrease 
rapidly  enough  to  refract  low  incident  angle  radio  waves  downward 
relative  to  the  earth  and  create  ducting  conditions.  If  the  decrease 
in  refractive  index  is  great  enough  a  surface-based  duct  occurs  which 
will  trap  surface-originated  radar  signals  and  give  greatly  extended 
ranges  for  detection  of  near  surface  targets. 
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Figure  la.  Radar  wave  path  under  "standard"  atmospheric  conditions. 
Note  path  curves  downward  but  at  a  rate  less  than  the  earth's  curvature, 
Beyond- the-horizon  target  detection  is  not  possible. 


Figure  lb.  Radar  wave  path  under  ducting  conditions.  Path  curves  down- 
ward at  a  rate  exceeding  the  earth's  curvature  resulting  in  beyond- the- 
horizon  target  detection. 
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The  air  adjacent  to  the  ocean  is  saturated  and  the  relative 
humidity  is  near  100%,  decreasing  rapidly  in  the  first  few  meters  to 
ambient  values  which  depend  on  varying  meteorological  conditions .  The 
decrease  of  humidity  causes  the  modified  refractive  index  (M)  to 
decrease  with  height  initially.  At  greater  heights,  the  humidity 
gradient  decreases  which  will  cause  M  to  reach  a  minimum  and  thereafter 
increase  with  height  (Fig.  2).  The  height  at  which  M  reaches  a 
minimum  value  is  called  the  evaporation  duct  height  (Z*)  and  is  a 
measure  of  the  strength  of  the  evaporation  duct. 
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Figure  2.  Relative  humidity  and  modified  refractivity  M  versus 
altitude  for  an  evaporation  duct. 
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The  probability  of  occurrence  of  duct  heights  large  enough  to 
cause  beyond- the-horizon  detection  capability  for  a  particular  radar 
varies  according  to  geographic  location,  season  and  the  time  of  day. 
The  strongest  ducts  are  generally  found  at  tropical  latitudes,  during 
the  summer  season,  and  during  daylight  hours. 
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III.  EVAPORATION  DUCT  MODEL 

Flux-profile  relationships  in  the  atmospheric  surface  layer 

* 

(Z<50m)  above  the  ocean  surface  have  been  studied  extensively  and  are 
well  understood.  Bulk  methods  can  be  used  to  determine  surface  flux 
scaling  parameters  from  routine  shipboard  measurements.  These  scaling 
parameters  can  be  applied  to  refractive  equations  to  determine  Z*. 
The  method  requires  only  mean  meteorological  conditions  as  input  para- 
meters and  is  well  suited  to  assessments  using  routine  shipboard  obser- 
vations . 

A.  EM  PROPAGATION  CONSIDERATIONS 

The  refractivity  N  (in  N  units)  of  the  atmosphere  for  EM  waves  is 
given  by  Battan  [1959] 

N  =  77.6  P/T  +  3.  73x10 5  e/T2  (1) 

where  P  is  the  atmospheric  pressure,  T  is  the  absolute  temperature 
and  e  is  the  partial  pressure  of  water  vapor.  An  expression  for  the 
vertical  gradient  in  N  is  given  by 

dN/dZ  =  CI  5P/6Z  +  C2  5q/5Z  +  C3  5T/6Z  (2) 

where  q  is  the  water  vapor  mixing  ratio  (q=.625e/P  in  gm/kg,  P  in  mb, 
T  in  °K,  and  Cl=.3,  C2=7.2  and  C3=-1.3).  A  negative  gradient  in  refrac- 
tivity at  height  Z  causes  a  horizontally  propagating  EM  wave  to  be  bent 
toward  the  surface  and  if  the  gradient  is  large  enough,  the  amount  of 
bending  exceeds  the  curvature  of  the  earth  and  the  wave  becomes  ducted. 
Since  magnitudes  of  specific  humidity  gradients  decrease  with  increas- 
ing height  in  the  near- surface  layer,  the  gradient  is  no  longer  strong 
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enough  for  ducting  at  some  height  Z*.  This  height  is  termed  the  evapo- 
rative duct  height  and  is  defined  by  the  critical  gradient  necessary 
for  EM  trapping. 

diN/dZ  =  -.157  (m_1) 
using  Equation  (2)  and  noting  dP/dZ=  -pq  (the  hydrostatic  balance  for 
the  atmosphere),  Z*  is  then  defined  by  the  following  equation. 

-.157  =  -.032  +  C2  5q/5Z  +  C3  6T/6Z  (3) 

Since  5q/6Z  and  5T/6Z  are  height  dependent  (approximately  1/Z)  there 
exists  a  value  of  Z=Z*  such  that  Equation  (3)  is  satisfied.  The  bulk 
method  allows  one  to  relate  5q/6Z  and  6T/6Z,  which  are  needed  to  deter- 
mine Z* ,   to  the  surface  layer  fluxes  and,  hence,  to  mean  parameters. 

B.  BULK  METHOD 

The  theoretical  framework  for  determining  the  evaporation  duct 
height  comes  from  knowledge  of  surface  layer  scaling  parameters  using 
the  Monin-Obukhov  similarity  (MOS)  theory.  The  theory  relates  the  sur- 
face layer  profiles  of  temperature,  wind  speed,  water  vapor  and  turbu- 
lence to  the  surface  fluxes  of  momentum,  sensible  and  latent  heat. 
Vertical  gradients  of  wind  speed  (U) ,  potential  temperature  (T)  or  water 
vapor  mixing  ratio  (q)  are  given  by  the  following  general  expression 

6X/5Z  =  (XVgkZ)  <j>x(Z/L)  (4) 

Where  X  =  U,  T  or  q,  k  =  Von  Karman's  constant  (0.35),  a  is  a  dif- 
fusivity  constant,  X*  is  the  scaling  parameter  and  (f>(Z/L)  is  the  MOS 
stability  correction.  The  surface  fluxes  define  the  scaling  parameters 
within  the  flux-profile  relationships.  Scaling  parameters  for 
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velocity,  potential  temperature  and  specific  humidity  are  related  to 
the  normalized  fluxes  by: 


U*  =  C-u'w')2    (m/s) 


T*  =  (-T'w')/U*  (°K) 


q*  =  (-q'w')/U*   (g/Kg) 


where  u'w',  T'w'  and  q'w'  are  the  mean  turbulent  fluxes  of  momentum,  heat 
and  moisture  respectively.  A  scaling  length,  the  Monin-Obukhov  length 
(L) ,  is  used  to  account  for  stability  influences .  Stability  magnitudes 
are  based  on  the  height  scaled  by  L,  i.e.  Z/L.  Negative  and  positive 
Z/L  values  corresponded  to  unstable  and  stable  conditions.  L  is  defined 
by  the  scaling  parameters  as 

L  =  TU*2  /*(Kg(T*  +  .61q*T))  (5) 

The  bulk  method  is  based  on  integral  forms  of  Equation  (4) .  Integrating 
Equation  (4)  from  Z  =  Z  ^  to  some  reference  height,  Z,  an  expression 
for  the  value  of  a  parameter  at  the  measurement  height  is  obtained 
[Businger,  1973]  • 

X(Z)  =  X(ZqX)  +  X*/aK(ln  z/ZqX  -  *(Z/L))  (6) 

where  KZ/L)  is  the  profile  stability  function  and  subscript 

X  =  U,  T,  or  q.   (Note:  a  and  ip  depend  on  X) 

Utilizing  Equation  (6) ,  a  method  of  estimating  surfaces  fluxes  in- 
volves the  scaling  parameters  and  roughness  lengths,  Z  „.  Appropriate 
values  of  temperature,  relative  humidity  and  wind  speed  are  obtained 
at  both  the  sea  surface  and  at  some  reference  height,  Z.  Equation  (6) 
can  be  used  to  determine  scaling  as  a  function  of  the  differences. 
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U*  =  (kU)  /  (In  (Z/Zq)  -  ih  (Z/L)) 

T*  =  (T-Ts)  cue/  (ln(Z/ZoT)  -<|;2  (Z/L)) 

q*  =  (q-qs)  ok/   (ln(Z/ZQT)   ^2  (Z/L)) 
a  was  set  equal  to  1.55  for  T  and  q  from  Businger  et  al  [1971]  .  The 
quantities  Z  and  Z  ~  are  roughness  lengths  for  velocity  and  tempera- 
ture profiles.  It  is  convenient  to  rewrite  these  equations  in  the 
generalized  drag  coefficient  form 

X*  =  Cx%  (X(Z)  -  X(0))  (7) 

The  neutral  stability  drag  coefficient  can  be  defined  as 

Cj^  =  aK/ln(Z/ZoX)  (8) 

then  the  drag  coefficient  becomes 

CX  =  °XN  '   (1  "  ((CXN^  /ctK)  ^    (Z/D))2  (9) 

In  this  form,  Z  ^  is  related  to  the  neutral  drag  coefficient  by 

Cjq*  =  ok/  ln(Z/ZoX)  (10) 

where  roughness  length 

2 

ZoX  =  Z  exp  (aK  /CXN  )  (11) 

Equation  (5)  can  be  rewritten  to  give  the  NDS  stability  parameter, 

Z/L,  as 

(1  -  (C/2)  /k   C¥i(Z/L)))* 

™ (12) 

Z/L  =  (Z/L), 


0  (1  -  (cJ2)  /cue  (4>2  (Z/L))) 


^TN 


where 


u 


<gZ     TN 


C™2    ((T(Z)-T(O))  +  .18(q(Z)-q(0))) 

(13) 


(Z/L)o=   T 


CUN  U(Z) 
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Note  that  (Z/L)  is  an  initial  estimate  of  stability  in  terms  of 
differences  and  neutral  drag  coefficient.  Final  values  of  Z/L  and 
scaling  parameters  can  be  found  from  a  simple  iteration  process . 
Stability  correction  functions  used  were: 
Velocity  profile: 

ipi(Z/I0  =  2  In  ((1  +  X)/2)  +  In  ((1  +  X  )/2) 

-2  ARCTAN   (X)    +  IT/2  Z/L<0 

X  =   CI  -15  Z/L)* 

ihCZ/L)  =  -4.7(Z/L)  Z/L>0 

Temperature  profile: 

ty2 (Z/L)  =  2  In  ((1  +  X)  /2)  Z/L<0 

X  =  (1  -9  Z/L)*! 
^2 (Z/L)  =  -6.5 (Z/L)  Z/L>0 

The  velocity  drag  coefficient,  C,^,  in  a  near- neutral  atmosphere 
was  calculated  from  the  wind- speed  dependent  formulation  given  by 
Kondo  [1975  ] .  A  typical  value  of  C^  at  Z=10  m  was  1.3x10" 3  m/s 
which  yielded  Z  =6x10  '   m.  Wind  speeds  and  drag  coefficients  are 
given  in  Table  I. 
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TABLE  I 
C^  VERSUS  WIND  SPEED  AT  10m  [  KONDO,  1975  ] 


UCm/s)  CDN  x  103 


.3-2.2  1.08xU"-15 

2.2-5.0  .77  +  .086  x  U 

5.0-8.0  .87  +  .067  x  U 

8.0-25.0  1.2  +  .025  x  U 

The  temperature  and  moisture  drag  coefficients,  G™,  and  C  N,  at  Z=10m 

were  estimated  as  1.3x10  3 [Davidson  et  al,  1978].   For  aT=1.35 

[Businger,  et  al,  1971  ],  a  roughness  length  of  Z  T=2.0xl0  5  m,  was 

obtained,  which  is  assumed  to  be  independent  of  wind  speed. 

C.  Z*  EQUATION  IN  MOS  FORM 

In  the  surface  layer,  the  gradient  of  N  can  be  written  in  terms  of 
T*  and  q*  as 

dN/dZ  =  -.032  =  <j>(Z/L)/(aicZ)  (7.2q*  -  1.3T*)  (14) 

which  arises  from  Equation  (3).  The  duct  height,  Z*,  is  obtained  from 

Equation  (14)  by  setting  dN/dZ  =  -.157,  the  critical  gradient  for  ducting, 

-(7.2q*  -  1.3T*) 

Z*  =    <P(Z*/L)  (15) 

qckC.125) 

Given  q*,  T*  and  L,  Equation  (15)  was  solved  iteratively  to  obtain  Z*. 

No  ducting  occurs  unless 

7.2q*  -1.3T*<  0 
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IV.  DATA  ACQUISITION 

The  data  used  in  this  analysis  were  unique  in  that  they  represented 
accurate  measurements  of  the  atmospheric  surface  boundary  layer  on  a 
horizontal  scale  of  1-200  km.  All  data  were  gathered  from  scientific 
research  vessels  whose  primary  cruise  objectives  were  to  observe,  dis- 
tinguish and  quantify  fluxes  and  mixing  in  atmospheric  and  oceanic 
boundary  layers.  The  following  are  descriptions  of  the  measurement 
platforms,  their  locations,  the  data  accuracy  and  the  intercomparison 
of  data  between  ships. 

A.  PLATFORMS  AND  LOCATION 

Data  for  this  analysis  were  taken  during  the  JASIN  experiment  in  the 

summer  of  1978.  Data  from  four  ships  involved  in  this  experiment  were 
selected  on  the  basis  of  two  objectives: 

1.  Examine  as  much  of  the  period  as  possible. 

2.  Examine  variation  over  nearly  fixed  spatial  separation. 

The  four  ships  selected  were  the  METEOR,  the  JOHN  MURRAY,  the  GARDLINE 
ENDURER  and  the  HECLA,  Figure  3  illustrates  the  profile  and  sensor  loca- 
tions for  each  ship.  Hourly  observations  were  made  of  air  temperature 
(T) ,  wet  bulb  temperature  (Tw) ,  sea  (bucket)  surface  temperature  (Ts) 
and  wind  speed  (U) .  The  triangle  configuration  of  the  ship  locations, 
depicted  in  Figure  4,  provided  quite  good  areal  coverage.  The  ships 
were  positioned  approximately  200  km  apart  at  each  corner  of  the 
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Figure  3.  Ship  profiles  and  sensor  locations:  (A)-  U, 
(B)  -  T,Tw  and  (C)  -  Ts. 
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Figure  4.  The  JASIN  area  in  relation  to  the  British  Isles 
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triangle  and  one  ship  was  located  at  the  center.  Nominal  positions 


were: 


MA  -  59.00  N  12.5  W 

MB  -  60.25  N  10.5  W 

MC  -  60.25  N  14.5  W 

CENTER  -  59.77  N  12.5  W 


B.  DATA  COLLECTION 

The  subject  observations  were  recorded  on  log  sheets  and  plots  of 
three-hourly  values  were  prepared  by  JASIN  participants.  The  plots 
from  the  GARDLINE  ENDURER,  JOHN  MURRAY  and  HECLA,  were  the  primary 
source  of  data  used  in  this  study.  Air  temperature,  wet  bulb  tempera- 
ture, sea  surface  temperature  were  accurate  to  0.2  °C,  pressure  to 
1  mb  and  wind  speed  to  1  kt.  Additionally,  data  used  from  the  METEOR 
were  available  in  computer  printout. 

C.  DATA  QUALITY  AND  INTERCOMPARISON 

During  the  period  15-22  July,  on  19  and  20  August,  and  on 
5  September,  the  four  ships  were  positioned  for  an  intercomparison  of 
meteorological  measurements.  Hourly  (half -hourly  during  formal  inter- 
comparisons)  values  of  wind  speed,  atmospheric  pressure,  dry  and  wet 
bulb  air  temperatures  and  sea  surface  temperature  were  compared  between 
vessels  approximately  one  kilometer  apart.  These  data  were  initially 
evaluated  for  coding  errors  and  inconsistencies  resulting  from  position 
changes.  Mean  and  standard  deviations  of  each  variable  were  calculated 
for  each  participating  vessel.  Time- series  of  inter-platform  sensor 
disagreements  were  produced  from  mathematical  differences  between 
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observed  variables.  Interplatform  means  and  standard  deviations  of 
these  differences  were  calculated  and,  assuming  the  disagreements  to  be 
normally  distributed,  the  901  confidence  interval  of  the  estimates  of 
differences  were  calculated. 

Results  derived  from  this  procedure  were  time  records  of  inter- 
platform sensor  differences.  Constant  disagreement  throughout  JASIN 
implied  the  sensor  performed  in  a  consistent  manner  and  it  remained 
only  to  explain  these  errors  in  terms  of  known  disparities ,  such  as 
non-uniform  sensor  heights  or  instrument  bias.  A  variation  of  the 
interplatform  corrections  required  further  analyses  either  in  terms  of 
atmospheric  stability  or  in  terms  of  trends  of  sensor  drifts.  The 
latter  were  established  by  regression  techniques.  Inconsistencies 
unexplained  or  uncorrected  by  these  considerations  remained  as  errors 
in  the  four -ship  system. 

The  METEOR  data  were  selected  as  the  standard  to  which  data  from  the 
other  vessels  were  adjusted.  However,  the  METEOR'S  sea  surface  tempera- 
ture values  were  corrected  with  data  from  the  other  three  ships.  The 
METEOR  was  chosen  because  it  was  involved  in  more  intercomparisdns , 
because  it  had  a  preferred  method  for  station-keeping,  and  because  of 
its  proximity  to  other  JASIN  ships  and  buoys.  The  corrections  accord- 
ing to  measured  parameters  and  individual  ships  are  summarized  below 

[Macklin  and  Guymer,  1980]  . 
1.  Wind  speed  (m/s) 

METEOR  =  METEOR 

HECLA  =  0.89  (HECLA)  +  0.3 

MURRAY  =  MURRAY 

ENDURER  =  ENDURER 
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2.  Pressure  (mb) 

METEOR  =  METEOR 

HECLA  =  HECLA  +  .08  To  0600,  July  28 

=  HECLA  +  .40  From  0900,  July  28 

MURRAY  =  1.00640  (MURRAY)  -  6.23 

ENDURER  =  ENDURER  +  .38  To  1700,  August  19 

=  ENDURER  +  .38  After  1700,  August  19   (P  less  than  1000) 

=  ENDURER  +  .48  After  1700,  August  19   (P  greater  than  1000) 

3.  Dry  Bulb  Temperature  (°C) 

METEOR  =  METEOR 
HECLA  =  HECLA  -  0.3 
MURRAY  =  MURRAY  -  0.5 
ENDURER  =  ENDURER  -  0.5 

4.  Wet  Bulb  Temperature  (°C) 

METEOR  =  METEOR 
HECLA  =  HECLA  -  0.3 
MURRAY  =  MURRAY  -  0.3 
ENDURER  =  ENDURER  -  0.1 

5.  Sea  Bucket  Temperature  (°C) 

METEOR  =  METEOR  +  0.3  To  1300,  July  27 

=  METEOR      1300,  July  27  to  0100,  August  3 
=  METEOR  -  0.4  0100,  August  3  to  1200,  August  9 
=  METEOR      From  1200,  August  9 

HECLA  =  HECLA 

MURRAY  =  MURRAY 

ENDURER  =  ENDURER 
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V.  ANALYSIS  PROCEDURES 

Data  analysis  procedures  were  designed  to  determine  1)  the  repre- 
sentative Z*  values  for  the  summer  regime,  2)  meteorological  factors 
which  most  affected  Z*,  3)  the  frequency  of  occurrence  and  duration  of 
ducts,  and  4)  the  extent  of  horizontal  homogeneity  of  the  duct.  To 
achieve  this,  the  data  were  first  corrected,  then  used  as  input  to 
bulk  formulae  to  compute  Z*. 

A  data  file  was  constructed  of  the  most  representative  Z*  values 
and  corrected  Ts-T,  U  and  RH  at  three-hour  intervals.  Winds  were 
reduced  from  anemometer  height  to  temperature  sensor  heights  using  the 
stability- corrected  logarithmic  wind  profile  [  Businger  et  al,  1971  ]. 
The  periods  examined  and  sensor  heights  are  summarized  in  Table  II. 
(Dates  are  in  year -month- day- hour  format) 

TABLE  II 
SHIP  DATA  PERIODS  AND  MEASUREMENT  HEIGHTS 
SHIP        STARTING     ENDING    ANEMOMETER  HGT  (m)   Z(m) 


MURRAY 

78071312 

78072103 

11 

6 

METEOR 

78072106 

78080912 

23 

11 

MURRAY 

78080915 

78081012 

11 

6 

METEOR 

78081709 

78090506 

23 

11 

ENDURER 

78090509 

78090606 

15 

4 

Statistical  and  time- series  analyses  of  Z*  were  performed  along 
with  subjective  analyses  of  influencing  meteorological  conditions.  The 
latter  were  performed  with  individual  observations,  synoptic  weather 
maps,  and  satellite  imagery  (visual  and  infrared).  These  descriptions 
were  correlated  to  the  time- series  of  Z*  values.  Joint  probability 
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density  distribution  and  conditional  mean  function  results  were  chosen 
to  examine  correlations  of  at  least  two  weather  parameters  to  Z* 
values.  Finally,  comparisons  were  made  of  Z*  and  all  atmospheric 
measurements  to  assess  the  horizontal  homogeneity  of  the  duct  for 
periods  when  ships  were  spatially  separated  and  at  fixed  locations. 

VI.  PRESENTATION  AND  INTERPRETATION  OF  RESULTS 

A.  GENERAL  ANALYSIS  OF  SUMMER  PERIOD  (13  July  -  6  September) 

Figure  5  is  a  time- series  comparison  of  evaporation  duct  height 
(Z*),  air- sea  temperature  difference  (Ts-T) ,  wind  speed  (U) ,  and 
relative  humidity  (RH)  for  the  period  of  study,  13  July  to  6  September, 
with  the  exception  of  the  interval  between  1200  GMT,  10  August  and  0900 
GMT,  17  August.  This  time  series  was  formulated  from  at  least  one  ship 
located  in  the  JASIN  array  and  provided  a  near -continuous  analysis  of 
the  evaporation  duct  height  (Z*)  at  three-hour  intervals.  General 
conditions  during  the  period  will  be  described  in  the  following  sections. 
1.  Synoptic  Conditions  in  JASIN  Area 

Daily  weather  maps  for  the  Northeast  Atlantic  are  shown  in 
Appendix  A.  Julian  day  and  date  (day  of  the  month)  are  shown  on  each 
map  (see  also  Table  A-I).  Studies  of  the  daily  surface  pressure  charts 
and  hourly  observations  [  Guymer  and  Taylor,  1978  ]  showed  that  there 
were  essentially  four  synoptic  situations  affecting  the  JASIN  area 
during  the  entire  period. 
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a.  Moist  Anticyclonic  (A) 

The  periods  of  14-18  July  and  22-30  August  were  affected  by 
moist  westerlies  or  northwesterlies  on  the  northern  flank  of  an  anti- 
cyclone to  the  west  of  Ireland.  Weak  frontal  waves  were  embedded  in 
this  flow  and  occasional  incursions  of  drier  air  took  place.  Winds 
varied  between  3  and  12  m/s. 

b.  Mobile  Westerly  Conditions  (B) 

The  periods  of  19-21  July  and  30  August  -  2  September  were 
characterized  by  depressions  moving  through  the  area  bringing  cold  and 
dry  air.  The  disturbances  were  very  weak. 

c.  Cyclonic  Southwesterly  (C) 

During  the  periods  22-29  July;  11-21  August  and  2-5  Sep- 
tember the  depressions  became  slow  moving  to  the  south  of  Iceland  and 
fronts  or  troughs  from  the  southwest  affected  the  region.  Generally, 
the  fronts  were  occluded  and  slow  moving  by  the  time  they  reached  the 
area,  with  only  moderate  winds.  However,  winds  which  were  the  strong- 
est experienced  during  the  entire  JASIN  period  reached  speeds  in  excess 
of  15  m/s  on  22/26  July,  18/20  August  and  5  September.  The  air  was 
significantly  warmer  than  the  sea. 

d.  Blocked  (D) 

During  the  period  30  July  -  10  August  a  blocking  low  formed 
over  northwest  France  and  moved  very  slowly  northeast;  to  the  west  a 
ridge  persisted.  Low  level  northeaster lies  were  colder  than  the  sea 
surface  but  convection  was  inhibited  by  warm  air  aloft  associated  with 
weak,  slow  moving  occlusions  which  also  produced  occasional  drizzle. 
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Towards  the  end  of  the  period,  the  ridge  intensified  and  moved  eastward, 
e.  Special  Events 

Some  interesting  atmospheric  occurrences  were  noted  as 
events  (refer  to  events  marked  on  Figure  5)  and  are  listed  below. 

1.  Relatively  cold,  dry  air  (14,  20  July;  overnight  1-2  August; 
7-9,  29,  31  August). 

2.  Significant  moisture  flux  in  near-neutral  conditions  (30  July; 
23-25  August;  1  September). 

3.  A  warm  frontal  wave  on  25  August  and  a  cold  front  preceded  by 
convective  instability  on  28  August.  Surface  fluxes  were  very  small. 

4.  Winds  in  excess  of  15  m/s  (22  July;  18,  19  August). 

5.  Incursions  of  warm,  dry  air  (14,  18,  25,  28,  29  July). 

The  following  are  descriptions  of  events  which  caused  warm,  dry  air 
(not  a  normal  occurrence  over  the  open  ocean)  to  occur. 

(1)  Afternoon  Heating.  On  14  July,  both  the  MURRAY  and 
HECLA  reported  increasing- air  temperature  but  MURRAY  reported  three 
consecutive  wet  bulb  temperature  readings  markedly  lower  than  HECLA. 
On  18  July,  MURRAY  and  METEOR  had  similar  air  temperature  readings 
while  HECLA  was  1.4°C  lower.  All  three  measured  similar  values  for 
wet  bulb  temperature.  On  both  occasions,  winds  were  light  under 
prevailing  high  pressure. 

(2)  Subsidence  and  Continental  Trajectory.  On  25  July, 

a  weak  ridge  moved  over  the  area  as  a  low  west  of  Ireland  moved  north- 
east. Winds  backed  from  the  northwest  to  the  east -southeast.  Warm, dry 
air,  a  result  of  subsidence  from  earlier  ridging,  flowed  from  the 
British  Isles  and  influenced  the  area  for  six  hours  before  the  approaching 
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occlusion  became  dominant. 

(3)  Short  Wave  Troughs.  Throughout  28  and  29  July,  a  long 
wave  trough  was  oriented  north- south  from  Iceland.  Weak  short  wave 
impulses  on  both  days  provided  surges  of  dry  air,  each  lasting  nine  hours, 
while  generally  maintaining  warm  southerly  flow. 

2.  Sea  Surface  Temperature  (TS) 

Surface  water  characteristics  were  associated  with  the  North 
Atlantic  Drift,  an  extension  of  the  warm  waters  of  the  Gulf  Stream.  The 
mean  water  temperature  in  the  vicinity  of  the  JASIN  triangle  was  12.2°C. 
The  general  areal  distribution  of  the  surface  temperature  for  the 
sampled  periods  is  depicted  in  Figure  6  [Liu  §  Katsaros,  1978]. 
The  Northeast  corner  of  the  triangle  was  always  the  coldest  position 
(11.2°C).  The  Northwest  and  South  corners  of  the  triangle  were  always 
the  warmest  (between  12  and  13°C) .  A  large  tongue  of  cool  water 
(11.8°C)  from  the  north,  intruded  south  and  south- southeast  through 
the  center  of  the  triangle.  It  persisted  for  the  entire  period  result- 
ing in  minimal  temperature  variations  at  the  center  position  (about 
0.2°C).  The  most  pronounced  Ts  gradients  occurred  along  a  line  between 
the  NW  and  NE  stations. 

3.  Representative  Z*  Values 

A  statistical  analysis  of  observed  parameters  and  calculated 
values  from  714  independent  observations  was  performed  to  determine  a 
representative  Z*  value  (Table  III).  The  mean  evaporation  duct  height 
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was  4.2  meters  with  a  standard  deviation  of  3.9  meters.  This  height 
appears  low  when  compared  with  the  climatology  of  the  region.  Sweet 
[1979]  determined  an* evaporation  duct  climatology  at  Ocean  Station 
ALFA  (62N/33W),  1100  km  upstream  of  the  JASIN  array,  and  Ocean 
Station  INDIA  (58N/19W) ,  400  km  southwest  of  the  JASIN  array.  The 
following  statistics  are  a  comparison  of  the  JASIN  evaporation  duct, 
using  the  method  previously  described  and  climatology  which  used  bulk 
aerodynamic  formulae  and  a  method  originated  by  Jeske  [1971]  and  modi- 
fied by  Hitney  [1975  ]. 

ALFA  JASIN  INDIA 

4.2  7 

3.8  7 

1.9-5.7  4.1-9.8 
714  6758 

The  interquartile  range  (IQR)  represents  the  lower  and  upper  bounds  for 

the  middle  duct  heights  (251  of  the  ducts  have  heights  beloitf  the  lower 

number  and  75%  have  heights  below  the  upper  number)  and  can  be  viewed 

as  an  indicator  of  the  spread  or  variation  of  the  duct  height. 

TABLE  III 

STATISTICS  FOR  JASIN  PERIOD  13  JULY  -  6  SEP  1978 

(714  OBSERVATIONS) 

TS    T     TS-T   U    RH    Z/L    T*   U*    q*    Z* 
(°C)  (°C)   (°C)  (m/s)  (%) (°C)  (m/s)  (g/kg)  (m) 

MEAN     12.2  11.8   .5    6.6   90    -.1   -.02  .27   -.05  4.2 

SIGMA     .6   .8   .9    2.9    8    1.0    .9   .13    .03  2.9 


— 

MEAN  (m) 
.MEDIAN  (m) 
IQR  (m) 
#  OBS 

5.8 
5.7 
3.4-8.3 
7058 
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4.  Atmospheric  Factors  Affecting  Z* 

Since  Z*  is  a  computational  artifact  determined  by  four  meteor- 
ological parameters,  Z*  (Ts-T,  Rh,  U,  Z/L) ,  multivariate  analyses 
were  used  to  determine  the  degree  to  which  individual  variables  affected 
duct  strength.  Joint  probability  density  distributions  were  computed 
for  each  of  the  parametric  interrelationships  to  examine  the  pattern  of 
atmospheric  factors  as  they  are  related  to  the  occurrence  of  Z*  values. 

Based  on  probability  theory  [Batchelor,  1953]  and  utilizing  a 
technique  described  by  Holland  [ 1968] ,  the  distribution  of  Z*  was  ana- 
lyzed statistically  against  the  distribution  occurrence  of  a  combination 
of  meteorological  variable  such  as  Ts-T  versus  U,  (Figures  7-12). 
Individual  moments  and  distributions  for  each  variable  were  also  computed 
in  the  analysis.  The  following  data  were  tabulated: 

1.  Mean 

2 .  Variance 

3.  Standard  Deviation  (Sigma) 

4.  Skewness 

5.  Kurt os is 

6.  Dimensional  Class  Limits  by  Sigma/4  Intervals  From  -4.0  Sigma 
to  +4.0  Sigma 

7.  Density  of  Observations  for  Each  Occurrence  by  Sigma/4  Class 
Interval 

Individual  class  and  joint  class  intervals  were  computed  for 
each  observation  set.  Two  indices  were  computed  for  each  data  pair 
representing  its  position  in  an  array  (18  x  18)  of  joint  class  inter- 
vals (Sigma/2  x  Sigma/2)  centered  on  the  mean  values  of  the  two 
variables.  Class  intervals  were  determined  by  dividing  the  standard 
value  of  one  of  the  joint  variables  by  0.5.  Pairs  of  independent 
variables  were  then  used  to  compute  joint  probability  density  functions 
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Figure  7.  Conditional  mean  (plotted  numbers)  of  Z*  for 
joint -frequency  distribution  (delineated  with  solid  labelled 
isolines)  of  wind  speed,  (U)  and  air-sea  temperature  differ- 
ence, (Ts-T).  Integers  along  the  right  and  bottom  edges  are 
sums  of  occurrences  in  respective  rows  or  columns.  The  total 
sum  is  714  but  may  be  less  if  the  compared  variables  are 
truncated  when  out  of  range. 
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(JDF) .  Similarly,  a  conditional  mean  function  (OIF) ,  defined  as  the 
average  values  of  Z*  for  all  observation  times  when  the  two  independent 
variables  had  values  within  the  joint  intervals,  was  defined  for  each 
joint  class  interval.  A  subjective  analysis  of  the  JDF's  and  the  OIF's 
was  performed  to  distinguish  relationships  as  being  from  computational 
artifacts  or  from  accidental  occurrences. 

a.  Air-Sea  Temperature  Difference  (Ts-T) :  (Figures  7,  8  and  9) 
The  air  temperature  averaged  0.5°C  colder  than  the  sea  sur- 
face temperature.  Z*(Ts)  calculations  were  largely  influenced  by  the 
observing  ship's  location  in  the  JASIN  area  (Fig.  6).  Air  temperature 
variations  were  primarily  influenced  by  migratory  weather  systems.  The 
most  notable  positive  deviation  of  Z*  occurred  when  Ts-T  was  less  than 
zero  (-.5°C),  with  low  relative  humidity  (RH<85%)  and  stable  atmospher- 
ic conditions  (Z/L>0).  However,  under  these  conditions  duct  strength 
was  extremely  sensitive  to  humidity  and  stability  changes.  If  RH  in- 
creased Sigma (RH)  or  Z/L  decreased  Sigma  (Z/L)  then  the  evaporation  duct 
virtually  vanished.  Generally,  Ts-T  had  little  effect  on  Z*  when 
RH>90%.  As  Ts-T  approached  zero  or  became  positive,  more  positive  Z* 
deviations  occurred  for  a  wider  range  of  RH  e.g.  (50-90%).  The  poorest 
ducting  conditions  occurred  when  T  exceeded  Ts  more  than  0.5°C  and  RH 
was  in  the  high  901 's.  In  near-neutral  conditions,  as  Ts-T  increased 
Z*   tended  to  increase  fractionally.  Generally  then,  Ts-T  was  more  of  a 
determinant  than  stability,  except  in  the  special  situation  of  warm, 
moist  air  over  cooler  water. 
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b.  Relative  Humidity:  (Figures  9,  10  and  11) 

Average  humidity  was  90%  and  ranged  from  49%  to  100%. 
Significant  changes  in  humidity  were  directly  due  to  the  passage  of 
large-scale  weather  systems.  As  relative  humidity  dropped  to  80%, 
marked  increases  in  Z*  occurred.  The  largest  Z*  values  occurred  when 
dry  (RH<80%)  and  warm  air  overlay  cool  water  (Ts-T  <-.5°C).  For  the 
relative  humidity  range  (50-90%),  larger  than  average  Z*  values 
occurred  when  the  air  was  1.0°C  cooler  than  the  sea  surface.  Lowest 
ducts  occurred  in  warm,  moist  meteorological  conditions  (RH: 85-100%, 
Ts-T:-2  to  0). 

c.  Wind  Speed  (U) :   (Figures  7,  10  and  12) 

Wind  was  most  influential  on  Z*  values  in  low  wind  situa- 
tions. In  light  winds,  changes  in  atmospheric  stability  (Z/L)  drasti- 
cally changed  the  duct.  As  seen  in  Figure  12,  a  change  in  Z/L  by 
Sigma  (Z/L)  changed  Z*  from  3  to  40  meters  for  winds  near  5  m/s.  In 
moist  air  (RH>90%) ,  Z*  was  low  regardless  of  wind  speed.  In  dryer  air 
(RH<85%)  wind  caused  increased  Z*  as  U  increased  from  5  to  12  m/s. 

d.  Stability  (Z/L):   (Figures  8,  11  and  12) 

Average  stability  was  near-neutral.  The  greatest  deviations 
of  Z*  occurred  when  Ts-T<0  (warm  air)  and  light  wind  (U<5m/s) .  Other- 
wise changes  in  Z/L,  hence  changes  in  Z*,  were  small  as  winds  increased 
or  colder  air  was  advected  over  relatively  warmer  water  (Ts-T>0). 

Figure  13  illustrates  the  interrelationships  that  the 
combination  of  selected  parameters  had  on  Z*.  The  ranges  of  Z*   are 
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depicted  for  the  primary  and  secondary  extremes  encountered  in  the 
joint  probability  density  distribution.  With  the  exception  of  the 
dependent  parametric  interrelationship  of  stability  (Z/L)  and  air- sea 
temperature  (Ts-T) ,  all  parameters  caused  large  values  of  Z*.  The 
primary  extreme  range  was  a  reflection  of  the  special  situation  in 
which  large  values  of  Z*  were  generated  due  to  atmospheric  conditions 
where  winds  were  light  and  warm,  dry,  stable  air  overlay  relatively 
cooler  waters.  A  more  significant  indication  of  atmospheric  influence 
was  the  secondary  range.  This  range  reflected  more  accurately  the 
effects  due  to  characteristic  changes  in  atmospheric  variables  with 
synoptic  conditions  of  a  summer  regime. 

The  most  influential  meteorological  factors  were  air-sea 
temperature  and  relative  humidity,  followed  closely  by  the  interrela- 
tionship of  wind  speed  to  air- sea  temperature.  Because  the  secondary 
range  of  wind  speed  to  relative  humidity  was  small,  this  indicated 
that  large  positive  Z*  deviations  were  influenced  most  significantly 
by  air-sea  temperature.  Stability  only  became  a  significant  factor 
in  light  winds.  Table  IV  summarizes  the  relationship  of  duct  strength 
to  atmospheric  parameters. 
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TABLE  IV 
DUCT  STRENGTH  VS.  ATMOSPHERIC  PARANETERS 

RH  (I)        U  (m/s)   Ts-T  (°C)      Z/L 


STRONGEST 

<  80 

<  5 

<  -.5 

>  0 

STRONG 

<  85 

5-12 

>  0 

WEAK 

85-100 

-2  to  0 

>  0 

WEAK 

>  90 

>  0 

5.   Z*  OCCURRENCE  AND  DURATION 

Histograms  of  duct  occurrence  and  duration  of  duct  heights  are 
shown  in  Figure  14.  This  figure  indicates  the  occurrence  of  strong 
evaporation  ducts  during  the  summer  is  small  and  short-lived.  Duct 
heights  of  6m  or  less  occurred  78%  of  the  time.  Rarely  did  significant 
ducts  occur  for  more  than  several  hours.  Figure  5  shows  the  changing 
weather  conditions  throughout  the  summer  and  the  effect  on  duct  height. 
Rarely  did  one  weather  pattern  dominate  or  persist.  Consequently,  the 
occurrence  of  strong  ducting  was  irregular.  It  was  most  closely  asso- 
ciated with  atmospheric  conditions  that  caused  either  cold,  dry  air  or 
warm,  dry  air  to  move  over  the  area.  Incursions  of  relatively  cold, 
dry  air,  usually  from  the  northwest  (Figures  21- 23,  28-33  and  40-45), 
caused  rapid  formation  of  evaporation  ducts  (9 -12m)  which  rarely  lasted 
longer  than  9-12  hours.  Dry  air  did  not  last  long  over  the  area  before 
it  was  modified.  Incursions  of  warm  dry  air  were  evident  and  produced 
brief  (3-6  hrs)  but  strong  evaporation  ducts  (>8m) .  Conditions  less 
favorable  for  the  formation  of  ducts  were  associated  with  cyclonic 
southwesterly  flow  of  warm,  moist  air  (Figure  20) . 
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Figure  15  and  16  compare  the  distribution  frequency  of  Z*.  Of  the 
three  locations,  JASIN  indicates  a  higher  percentage  of  low  duct 
heights,  Ocean  Station  INDIA  shows  its  mode  (group  that  occurs  most 
frequently)  in  the  6-9  m  range,  the  highest  of  the  three.  Additionally, 
INDIA  had  a  greater  percentage  of  Z*  in  the  9- 18m  range.  Ocean  Station 
ALFA   located  off  the  east  coast  of  Greenland  and  in  an  oceanic  region 
of  confluence  where  the  cold  Greenland  current  meets  the  cool  Irminger 
current,  indicated  a  mode  of  3- 6m  and  a  generally  higher  percentage  of 
Z*  greater  than  JASIN  but  less  than  INDIA. 

This  difference  between  the  evaporation  duct  climatology  in  the 
three  areas  appears  to  have  been  due  to  (1)  the  proximity  of  ALFA  to 
the  Greenland  landmass  and  the  associated,  extremely  cold  water  and 
(2)  INDIA'S  location  in  the  deeper  and  warmer  water  of  the  North 
Atlantic  current.  Interpretation  of  visual  and  infrared  satellite 
imagery  indicates  (Figures  23-33)  the  flow  of  cold,  dry  air  off 
Greenland,  coupled  with  cold  sea  surface  temperatures,  causes  higher 
Z*  values  at  ALFA  due  to  the  large  difference  in  surface  humidity 
and  atmospheric  humidity.  As  the  dry  air  mass  moves  across  the  North 
Atlantic  atmospheric  humidity  changes  are  relatively  small  as  it 
approaches  both  the  JASIN  and  INDIA  locations.  The  key  difference 
between  the  two  locations  is  Ts.  The  warmer  water  at  INDIA  causes 
higher  surface  humidity,  thus  a  larger  difference  in  atmospheric  and 
surface  humidity  resulting  in  the  occurrence  of  stronger  ducts.  JASIN, 
because  the  air  mass  has  become  more  moist  since  leaving  ALFA,  but 
enters  an  area  of  relatively  warmer  water,  experiences  the  smallest 

44 


50 


40 


% 


30 


20 


10 


! 1 i 1 

9        12       15        18      21       OO 
Meters 


JASIN    AREA 

Figure  15.     Frequency  of  Z* 
occurrence   (13  Jul  -  6  Sep  1978) 


OCEAN  STATION  ALFA 


OCEAN  STATION  INDIA 


50 
40 
%  30 
20 
10 


50 
40 
%  30 
20 
10 


0  3  6  9  12  15  18  21  30  • 

Meters 

July 


T. 


0  3  6  9  12  15  18  21  30  - 

Meters 

July 


Figure  16.  Frequency  of  Z*  occurrence  (climatology)  at 
North  Atlantic  Ocean  Stations  nearest  JASIN, 


45 


Figure  17.  Ocean  Stations  in  the  North  Atlantic:  ALFA 
(62N,  33W)  and  INDIA  (58N,  19W) . 
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difference  in  atmospheric  and  surface  humidities,  therefore  the  weakest 
ducts  of  the  three  locations. 

B.  HORIZONTAL  NATURE  OF  THE  EVAPORATION  DUCT 

The  three  periods  (7-8,  23-25  August  and  30  August  -  1  September) 
examined  corresponded  to  times  when  ships  were  positioned  at  all  corners 
of  the  JASIN  triangle.  The  objectives  of  analysis  for  these  periods 
were  to: 

1.  Determine  the  horizontal  homogeneity  of  the  evaporation  duct. 

2.  Determine  which  meteorological  factors  enhanced  or  degraded  areal 
ducting. 

The  analysis  procedure  for  this  was  to: 

1.  Compare  computed  Z*  values  at  different  ship  positions  and  con- 
comitant meteorological  parameters  to  assess  similarity. 

2.  Analyze  synoptic  meteorological  conditions  utilizing  satellite 
imagery  (visual/ infrared) ,  ship  observations  and  synoptic  weather 
maps. 

3.  Analyze  the  temporal  variations  or  consistency  of  horizontally 
homogeneous  conditions. 

4.  Relate  the  above  fields  to  the  spatial  nature  of  the  evaporation 
duct. 

Area  homogeneity  was  determined  by  comparing  Z*,  Ts-T,  RH  and  U 

for  ships  positioned  at  each  corner  of  the  triangle.  Azimuthal  sectors 

were  defined  as  the  legs  of  the  triangle  (NW  sector:  NW-S,  NE  sector: 

NE-S,  E  sector:  NE-NW) .  Using  the  South  (S)  position  as  the  reference, 

differences  were  determined  relative  to  Northwest  (NW)  and  Northeast 

(NE)  locations  (NW  and  NE  sectors,  respectively).  Similarly,  using  the 

Northwest  observer  as  reference,  differences  were  determined  relative 

to  the  Northeast  location  (East  sector) .  Differences  between  locations 

(ships)  relative  to  sectors  are  tabulated  in  Tables  V  through  VII. 
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1.  Z*  Range  in  Homogeneous  Atmosphere 

A  necessary  step  to  interpret  the  results  of  the  tabulated 
differences  is  the  performance  of  an  error  analysis  to  determine  ttfe 
possible  range  of  computed  Z*  which  could  occur  in  a  homogeneous 
atmospheric  situation  due  to  measurement  inaccuracies.  An  assessment 
was  required  to  ascertain  the  maximum  deviation  allowable  in  Z*  for  the 
three- ship  system  under  nearly  homogeneous  atmospheric  conditions. 
Additionally,  a  sensitivity  analysis  of  Z*  values  was  performed  for 
each  inherent  measurement  uncertainty. 

Previous  studies  on  bulk  measurements  [Davidson  et  al,  1978] 

have  concluded  that  representative  inaccuracies  in  bulk  parameter 

measurements  are: 

RH  =  1% 
Ts-T  =  0.5° 
U  =  10% 

Five  observation  times  were  chosen  to  perform  the  error  analysis 

These  corresponded  to  times  when  air  temperature,  relative  humidity  and 

wind  speed  were  nearly  similar  for  all  ships  stationed  on  the  triangle. 

Benchmark  Z*   values  were  computed  from  each  of  these  representative 

homogeneous  atmospheric  situations.  Atmospheric  conditions  at  1500  GMT 

August  23,  24  and  0600  GMT  25  August  were  similar  at  all  three  ship 

locations.  Mean  values  of  bulk  parameters  were  used  to  compute  a  Z* 

at  each  corner  of  the  triangle  assuming  no  measurement  inaccuracies. 

A  similar  analysis  was  performed  on  data  taken  at  0600  and  1500  GMT  31 

August.  The  following  gives  the  results  of  the  two  periods: 
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CASE  I 

Ts  =  12  °C        T  =  11.8  °C 

SHIP        Z(m)         Z/L 

ENDURER        4  -.0055 

MURRAY         6  -.0082 

METEOR        11  -.0139 

CASE  II 


RH  =  901 

U  =  8.1  m/s 

Z*(m) 

4.44 

4.32 

4.17 

Ts  =  12.5  °C 

T  = 

11.5°C 

SHIP 

Z(m) 

Z/L 

RH  =  75%    U  =  10  m/s 
Z*(m) 


ENDURER        4  -.015        10.4 

MURRAY         6  -.023        10.1 

METEOR        11         -.045        9.7 


The  Z*  values  computed  in  each  of  the  above  cases  represented 
reference  Z*  values  for  each  measurement  height,  Z  (e.g.  4,  6,  11m)  in 
a  spatially  homogeneous,  mixed  layer.  Measurement  uncertainties  listed 
below  were  introduced  into  the  Z*  calculations  and  the  maximum  error 
computed.  Percentage  values  for  both  cases  reflect  the  change  in  Z* 
due  only  to  the  measurement  errors  of  each  individual  parameter. 

CASE  I  CASE  II      .  MEASUREMENT  ERROR (rms) 

3%  10% 

81  3% 

2%  .5  °C 

The  cumulative  error  in  the  worst  case  situation,  where  measurement 

uncertainties  were  between  two  ships  in  a  horizontally  homogeneous 

atmosphere,  produced  the  following  difference  in  Z*  values: 
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u 

7% 

RH 

25% 

Ts 

14% 

Z*  Cm)  Z*  Diff  (Max-Min) 

4  2.7 

10  2.9 

Therefore,  it  was  concluded  that  if  the  atmosphere  were  homo- 
geneous and  the  only  uncertainties  were  in  measurements  of  T,  Ts,  RH 
and  U,  then  the  greatest  difference  in  Z*  that  any  two  ships  would 
have  experienced  was  2.9  m. 

2.  Case  Studies  of  Horizontal  Homogeneity 

Three  periods  were  examined,  in  view  of  the  error  analysis,  when 
ships  were  stationed  at  the  corners  of  the  JASIN  triangle  (approximately 
200  km  apart).  Simultaneous  meteorological  observations,  taken  at  three- 
hour  intervals,  and  their  respective  Z*   values  were  compared.  In  each 
case,  an  average  or  better  than  average  evaporation  duct  existed.  Exten- 
sive periods  of  below  average  ducting  were  not  considered. 

The  examination  included  consideration  of  the  general  synoptic 
conditions,  differences  in  local  meteorological  conditions  between  ships 
(and  the  factors  which  most  influenced  Z*) ,  the  extent  of  the  duct  over 
the  entire  triangle  (areal  homogeneity)  and  between  corners  (sector  homo- 
geneity) .  The  maximum  Z*  difference  (2.9m)  between  ships  (due  to 
measurement  errors  as  previously  described)  was  the  criterion  for  duct 
homogeneity.  Therefore,  if  the  maximum  difference  in  Z*  between  all 
three  ships  was  less  than  or  equal  to  2.9m,  areal  homogeneity  was 
assumed  to  exist.  Similarly,  if  the  criterion  was  valid  between  corners 
of  the  triangle,  sector  homogeneity  was  assumed  to  exist. 
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a.  Case  I 

(1)  Period.   7  Aug  (0000  GMT)  -  8  Aug  (2100  GMT) 

(2)  Observer  Positions  on  Triangle. 

Northwest  (NW)  -  ENDURER 
Northeast  (NE)  -  MURRAY 
South  (S)      -  METEOR 

(3)  Synoptic  Conditions  -  Blocked.  A  ridge  persisted  over 
the  area  as  it  intensified  and  moved  slowly  east  (Figure  19) .  Northerly 
surface  winds  were  moderate  in  the  eastern  region  and  lighter  in  the 
west.  Low  level  instability  prevailed  as  cold  air  was  advected  over 
warmer  water.  However,  convection  was  inhibited  by  warm  air  aloft  as 

an  occlusion  moved  slowly  over  the  area  (Figures  25-27). 

(4)  Results  of  Tabulated  Differences  -  (Table  V) .  Obser- 
vations at  NW  and  NE  were  more  humid  than  S.  Winds  were  stronger  at 
NE  than  NW  and  S.  Air-sea  temperature  differences  were  greatest  at  S. 
The  duct  heights  on  the  average  were  2m  higher  in  the  south  than  the 
two  northern  positions  (i.e.  6.9m  vs.  4.9m). 

(5)  Areal  Homogeneity.  Similar  evaporation  duct  strengths 
occurred  simultaneously  at  the  three  stations  30%  of  the  time  with  the 
most  sustained  conditions  lasting  12  hours.  Figure  18  illustrates  by 
sectors  the  percentage  of  homogeneous  ducts  and  the  longest  duration 
experienced.  Northerly  flow  coincided  with  similarity  between  the  two 
northern  positions  while  the  greatest  disparity  occurred  between 
southern  and  western  positions.  The  most  persistent  horizontal  homo- 
geneity occurred  between  NW  and  NE  where  humidities  and  air -sea 
temperature  differences  were  most  similar.  Poorest  conditions  existed 
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between  NW  and  S  due  primarily  to  the  dissimilarity  in  Ts-T  and  RH  and 
the  consequent  difference  in  air-sea  humidity  (Figures  44-47). 
b.  Case  II 

(1)  Period:  22  August  (1500  GMT)  -  25  August  (1500  GMT) 

(2)  Observer  Positions  on  Triangle: 

Northwest  (NW)  -  ENDURER 
Northeast  (NE)  -  HECLA 
South  (S)      -  METEOR 

(3)  Synoptic  Conditions  -  Moist  Anticyclonic .  An  advanc- 
ing cold  front  passed  over  the  area  the  afternoon  of  23  August  (Fig- 
ures 21,  28,  31).  Cold,  dry  northwesterly  air  immediately  followed  the 
frontal  passage  and  prevailed  through  24  August  (Figures  22,32,  33). 
Humidity  gradually  increased  from  the  west  over  the  area  by  25  August 
as  a  warm  front  approached  (Figures  23,  34,  35). 

(4)  Results  of  Tabulated  Differences  -(Table  VI).  Generally, 
the  greatest  differences  experienced  were  between  NW  and  NE  stations 

due  to  the  west  to  east  movement  of  the  air  mass  over  the  area.  Winds 
and  humidity  were  similar  at  all  observing  stations.  Air-sea  differ- 
ences were  least  at  the  NE  position  where  water  was  the  coldest,  hence 
weakest  ducts  were  generated.  Duct  heights  ranged  1.1-12. 2m  and 
averaged  5m  over  the  area. 

(5)  Areal  Homogeneity.  The  area  experienced  homogeneous 
ducting  44%  of  the  time  with  the  longest  duration  9  hours.  Homogeneity 
was  nearly  equal  for  all  sectors  with  the  northeast  experiencing 
slightly  stronger  ducting  and  lasting  for  24  hours  following  frontal 
passage  (Figure  18) .  The  greatest  Z*  difference  across  the  front  was 
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3.8m  i.e.  6.6m  (NW)  vs.  2.8m  (NE)  Horizontally  homogeneous  ducting 
lasted  for  nearly  15  hours  following  frontal  passage  before  another 
impulse  of  dry  air  entered  the  western  area.  It  took  15  hours  for  the 
dry  air  to  move  across  the  JASIN  area  and  strengthen  the  duct  in  the 
east  (Figures  48-55) . 
c.  Case  III 

(1)  Period.  30  August  (0000  GMT)  -  1  September  (2100  GMT) 

(2)  Observer  Positions  on  Triangle. 

NORTHWEST  (NW)  -  ENDURER 
NORTHEAST  (NE)  -  HECLA 
SOUTH  (S)      -  METEOR 

(3)  Synoptic  Conditions  -  Mobile  Westerly.  A  weak 
depression  moved  over  the  region  (Figure  19).  The  JASIN  area  was  in  a 
warm  sector  on  30  August  as  a  cold  front  passed  early  on  31  August.  A 
warm  front  approached  from  west  late  on  1  September  (Figures  40-43.2). 

(4)  Results  of  Tabulated  Difference  -  (Table  VII) .  The 
average  duct  height  was  6m.  The  most  intense  ducts  occurred  on  31 
August  after  the  cold  front  passed.  They  occurred  in  the  west  and 
south  where  sea  surface  temperatures  were  highest  and  air- sea  tempera- 
ture differences  were  greatest.  On  the  average,  relative  humidity 
was  lower  to  the  northeast  than  the  western  and  southern  region  which 
were  largely  under  the  influence  of  warm  sector  of  the  low. 

(5)  Areal  Homogeneity.  Due  to  the  high  mobility  of  the 
weather  systems,  areal  horizontal  homogeneous  ducting  occurred  only  29% 
of  the  time  with  a  maximum  duration  of  9  hours  during  very  low  ducts. 
Sector  homogeneity  was  poorest  between  the  northeast  and  northwest 
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positions  and  best  between  south  and  northeast  (Figure  18) .  The  maxi- 
mum difference  in  Z*  across  the  cold  front  was  9.5m  i.e.  11.7m  (NW) 
vs.  2.2m(S).  The  strongest  duct  heights  associated  with  any  weather 
system  were  generated  on  31  August  following  the  cold  front.  Though 
outside  the  limits  of  horizontal  homogeneity  criterion,  2.9m,  duct 
heights  of  8. 6-13. 4m  occurred  at  all  stations  for  12  hours  before 
weakening  in  the  east  (Figures  54-59). 

3.  Summary  of  Horizontal  Variation  of  Z* 

The  periods  of  the  three  cases  studied  were  chosen  because 
higher  than  normal  ducts  occurred  for  reasonably  long  periods.  The 
horizontal  nature  of  significant  ducting  in  the  region  could  be  assessed. 
Each  case  was  unique  in  that  strong  ducting  was  the  result  of  differing 
synoptic  meteorological  conditions. 

The  strongest  ducts  occurred  as  mobile  westerly  depressions 
moved  through  the  area;  however,  horizontal  homogeneity  was  poorest. 
Air- sea  temperature  differences  and  relative  humidity  were  the  most 
important  parameters  affecting  duct  strength.  Larger  Z*   values  were 
experienced  at  the  NW  and  S  locations  because  these  locations  had  higher 
values  of  Ts  than  the  NE  location.  Since  Ts  gradients  did  exist  between 
observation  positions,  sector  horizontal  homogeneity  appeared  to  be 
influenced  by  the  direction  from  which  the  air  mass  entered  the  area. 
The  NW  and  NE  positions,  which  experienced  the  largest  Ts  gradient,  had 
more  homogeneous  ducting  with  northerly  flow  than  westerly  flow.  The 
gradient  between  S  and  NE  was  small  and  resulted  in  consistently  homo- 
geneous conditions  in  all  three  cases. 
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Case  I  illustrated  the  effect  wind  plays  on  Z*  values.  Both 
NE   and-NW  had  high  relative  humidity,  but  NW  had  a  greater  air-sea 
temperature  difference,  resulting  in  a  larger  moisture  flux.  However 
the  stronger  winds  in  the  eastern  region  compensated  for  the  lower 
moisture  flux  producing  nearly  similar  duct  heights.  A  comparison  of 
the  S  and  NE  stations,  which  had  similar  winds  and  air-sea  temperature 
differences,  showed  significantly  different  duct  heights.  The  import- 
ant parameters  in  this  case  were  relative  humidity  and  sea  surface 
temperature  which  produced  larger  surface  moisture  fluxes. 

Case  II  exhibits  the  effects  due  to  incursions  of  cold,  dry  air 
with  moderate  to  fresh  winds.  Satellite  pictures  (Figures  50,  51) 
indicated  low  level  instability  under  windy  and  drying  conditions.  How- 
ever, significant  afternoon  heating  sharply  reversed  the  air-sea  tempera- 
ture difference  causing  the  ducts  to  strengthen.  Despite  the  afternoon 
variation,  this  case  had  the  most  consistent  horizontally  homogeneous 
ducting,  occurring  44%  of  the  time,  as  well  as  the  most  equal  sector 
homogeneity  of  the  three  cases.  The  strongest  ducts  were  toward  the 
northwest  where  the  air  was  the  driest.  The  strong  zonal  flow  in- 
hibited the  driest  from  reaching  the  southern  region  where  relative 
humidity  varied  little  and  air-sea  temperature  differences  were  small. 

Case  III  illustrates  the  effect  of  rapidly  moving  depressions 
on  the  formation  and  duration  of  ducts.  The  strongest  areal  ducts 
formed  with  the  passage  of  a  cold  front  but  lasted  only  a  short  period 
before  the  next  upstream  weather  system  weakened  the  duct  in  the 
western  region.  The  occurrence  of  ducts  in  the  NE  lagged  the 
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occurrence  of  similar  ducts  at  the  NW  position  by  six  hours. 

A  composite  of  the  frequency  of  occurrence  of  areal  and  sector 
homogeneous  ducts  and  maximum  duration  is  depicted  in  Figure  18.  It 
was  concluded  that  areal  homogeneity  was  a  probable  occurrence  only  35% 
of  the  time  for  the  three  periods  studied,  with  a  maximum  likely  dura- 
tion of  12  hours.  The  NE  sector  appeared  to  maintain  the  best  horizon- 
tal homogeneous  conditions.  The  sector  also  appeared  to  exhibit  good 
homogeneity  but  of  short  duration.  The  NW  sector  exhibited  the  least 
homogeneous  conditions  but  when  homogeneity  occurred,  it  lasted  up  to 
24  hours. 
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Figure  18.  Percentage  of  areal  and  sector  horizontal 
homogeneity  of  ducting  and  maximum  duration  for  three 
cases  and  the  average. 
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VII.  SUMMARY  AND  CONCLUSIONS 

Unique  temporal  and  spatial  descriptions  of  evaporation  duct 
features  in  the  Northeast  Atlantic  were  described  from  observational 
data  taken  from  four  ships  spaced  in  a  fixed  array.  Evaporation  duct 
heights  were  determined  from  accurately  measured  bulk  surface  data 
utilizing  Monin-Obukhov  scaling  in  flux-profile  relationships.  Discern- 
ible synoptic  features  from  surface  pressure  fields  and  DMSP  satellite 
pictures  were  air  mass  histories,  trajectories  and  horizontal  variations 
(frontal  zones) .  Synoptic  scale  spatial  variations  of  the  sea  surface 
temperature  were  also  available  and  were  as  important  in  the  spatial 
variation  as  the  atmosphere. 

Evaporation  duct  features  established  for  this  data  set  were: 

1.  The  mean  Z*  values  were  less  than  those  reported  in  available 
climatologies  based  on  weather  ships'  observations  (4.2m  versus 
7.0m). 

2.  Z*  values  at  or  above  8-10  meters  occur  frequently  in  conjunction 
with  outbreaks  of  cold,  dry  air  from  the  north  or  northwest  and 
have  durations  of  9  to  12  hours. 

3.  Air-sea  temperature  differences  and  relative  humidity  were  the 
important  surface  layer  parameters  in  establishing  tactically 
significant  values  and  the  former  was  the  most  important. 

4.  If  horizontal  homogeneity  were  defined  as  representatively  uniform 
Z*  values  over  200  km  distances,  such  conditions  occurred  only 
35%  of  the  times  when  significant  duct  heights  existed,  and  for 
durations  of  up  to  12  hours. 

Conclusions  relative  to  synoptic  scale  features  most  closely  asso- 
ciated with  these  near- surface  features  are: 

1.  The  primary  synoptic  influence  on  the  duct  height  is  the  air-mass 
trajectory  (due  to  its  role  in  determining  the  moisture  content) 
and  the  air- surface  temperature  difference.  North  and  northwest 
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trajectories  have  the  highest  associated  Z*  values. 

2.  The  primary  feature  influencing  horizontal  homogeneity  over  dis- 
tances of  200km  is  the  sea  surface  temperature,  and  its  influence 
depends  on  the  air  mass  trajectory. 

3.  Light  wind  conditions  modify  the  above  specified  synoptic  feature- 
surface  duct  feature  relationships.  The  modifications  occur  over 
smaller  temporal  and  spatial  scales  than  those  of  the  prevailing 
features . 

A  concluding  assessment  is  that  tactically  significant  features  of 
the  evaporation  duct  in  the  study  area  and  for  the  study  period  can  be 
established  from  synoptic  scale  descriptions.  Further  exploitation  of 
available  data  should  emphasize  satellite  imagery  which  could  delineate 
surface  temperature  and  air  mass  characteristics. 

A  recommendation  for  further  efforts  is  that  all  available  surface 
data  from  the  JASIN  project  be  used  in  order  to  extend  and  substantiate 
feature  descriptions  and  relationships  identified  in  this  study.  Since 
the  start  of  this  study,  additional  surface  layer  data  from  buoys  and 
other  participating  ships  have  been  compiled  by  the  organizing  institu- 
tions and  made  available  on  computer  tapes.  These  data  would  increase 
the  observations  used  by  two  orders  of  magnitude,  thereby  improving  the 
descriptions  of  both  the  temporal  and  spatial  features. 
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APPENDIX  A. 
TABLE  A- I 


Julian   (year)   Dav  and  Date 


Date 

Day 

Dace 

Day_ 

Date 

Day_ 

July  12 

193 

Aug.   1 

213 

Sept.   1 

244 

13 

194 

2 

214 

2 

245 

14 

195 

3 

215 

3 

246 

15 

196 

4 

216 

4 

247 

16 

197 

5 

217 

5 

248 

17 

198 

6 

218 

6 

249 

18 

199 

7 

219 

7 

250 

19 

200 

8 

220 

3 

251 

20 

201 

9 

221 

9 

252 

21 

202 

10 

222 

10 

253 

22 

203 

11 

223 

11 

254 

23 

204 

12 

224 

12 

255 

24 

205 

13 

225 

13 

256 

25 

206 

14 

226 

14 

257 

26 

207 

15 

227 

15 

258 

28 

209 

16 

228 

16 

259 

29 

210 

17 

229 

17 

260 

30 

211 

18 

230 

31 

212 

19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 

231 
232 
233 
234 
235 
236 
237 
238 
239 
240 
241 
242 
243 
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Figure  20,  500  mb  (Top)  and  SFC  Pressure  Maps  (Bottom) , 
1200  GMT  19  Aug  78. 
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Figure  21.  500  mb  (Top)  and  SFC  Pressure  (Bottom)  Maps,  1200  GMT, 


23  Aug  78, 
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Figure  22. 
24  Aug  78. 


500  mb  (Top)  and  SFC  Pressure  (Bottom)  Maps,  1200  GMT, 
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Figure  23. 
25  Aug  78. 


500  mb  (Top)  and  SFC  Pressure  (Bottom)  Maps,  1200  GMT, 
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APPENDIX  B 


Defense  Meteorological  Satellite  Program  (DMSP)  imagery  (visual 
and  infrared)  of  the  Northeast  Atlantic  and  high-resolution  imagery 
of  the  JASIN  area  corresponding  to  periods  of  significant  ducting. 
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Figure  24.  DMSP  visual  imagery,  25  July,  1212  GMT. 
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Figure  25.  DMSP  infrared  imagery,  7  August,  1032  GMT 
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Figure  26.  EMSP  visual  imagery,  8  August,  1143  GMT. 


Figure  27.  DNBP  infrared  imagery,  8  August,  1143  GMT, 
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Figure  28.  EMSP  visual  imagery,  22  August,  1110  (MI, 


Figure  29.  DMSP  infrared  imagery,  22  August,  1110  GMT, 
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Figure  30, 


EMSP  visual  imagery,  23  August,  1221  CMT. 


Figure  31.  EMSP  infrared  imagery,  23  August,  1221  GNfT. 
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Figure  32.  EMSP  visual  imagery,  24  August,  1138  GNfT. 


Figure  33.  EMSP  infrared  imagery,  24  August,  1138  GMT, 
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Figure 


34.  EMSP  visual  imagery,  25  August,  1054  GMT. 


Figure  35.  EMSP  infrared  imager)',  25  August,  1054  GMT, 
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Figure  36.  EMSP  visual  imagery,  28  August,  1233  GMT. 


Figure  37.  DMSP  infrared  imagery,  28  August,  1233  GMT, 
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Figure  38.  DMSP  visual  imagery,  29  August,  1150  GMT, 
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Figure  39.  EMSP  infrared  imagery,  29  August,  1150  QW , 
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Figure  40.  EMSP  visual  imagery,  30  August,  1107  GMT, 


Figure  41.  EMSP  infrared  imagery,  30  August,  1107  GMT. 


81 


':3sip^>^ 


Figure  42.  DMSP  visual  imagery,  31  August,  1218  GMT, 


Figure  43.  DMSP  infrared  imagery,  31  August,  1218  GMT, 
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Figure  43.1  EMSP  visual  imagery,  1  Sept,  1135  GMT. 


Figure  43.2  DMSP  IR  imagery,  1  Sept,  1135  GMT. 
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Figure  44.  High- resolution  DMSP  visual  imagery  and  surface  data. 
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High- resolution  DMSP  IR  imagery  and  surface  data, 
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Figure  46 


High- resolution  DMSP  visual  imagery  and  surface  data. 
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Figure  47.  High- resolution  DMSP  IR  imagery  and  surface  data. 
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Figure  48.  High- resolution  DMSP  visual  imagery  and  surface  data. 
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Figure  49.  High-resolution  MSP  IR  imagery  and  surface  data. 
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Figure  50.  High- resolution  DMSP  visual  imagery  and  surface  data, 
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Figure  51.  High -resolution  DMSP  IR  imagery  and  surface  data. 
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Figure 


52.  High-resolution  DMSP  visual  imagery  and  surface  data. 
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Ficoire  53.  Hi gh- resolution  DMSP  IR  imagery  and  surface  data. 
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Figure  54.  High -resolution  EMSP  visual  imagery  and  surface  data, 
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Figure  55.  High- resolution  DMSP  IR  imagery  and  surface  data. 
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Figure  56.  High-resolution  DMSP  visual  imagery  and  surface  data. 
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Figure  57.  High- resolution  DMSP  IR  imagery  and  surface  data. 
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Figure  58.  High- resolution  DMSP  visual  imagery  and  surface  data 
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Figure  59-  High-resolution  DMSP  IR  imagery  and  surface  data, 
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